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ABSTRACT 

We study dynamical capture binary neutron star mergers as may arise in dense stellar regions such 
as globular clusters. Using general-relativistic hydrodynamics, we find that these mergers can result 
in the prompt collapse to a black hole or in the formation of a hypermassive neutron star, depending 
not only on the neutron star equation of state but also on impact parameter. We also find that these 
mergers can produce accretion disks of up to a tenth of a solar mass and unbound ejected material of 
up to a few percent of a solar mass. We comment on the gravitational radiation and electromagnetic 
transients that these sources may produce. 

Subject headings: black hole physics — gamma-ray burst: general — gravitation — gravitational waves — 
stars: neutron 



1. INTRODUCTION 

Merging binary neutron stars (NSs) warrant de- 
tailed study because these systems promise to be 
rich sources of both gravitational and electromagnetic 
(EM) radiation, probing strong-field gravity and nu- 
clear density physics. NS-NS mergers are a pri- 
mary source targeted by g ravitational wave (GW ) 
detectors (such as LIGO; lAbramovici et al.l I1992D. 
They are also candidates for short gamma-ray burst 
(SGRB) progenitors and several other EM coun- 
terparts (iMetzger fc Bergerl I2Q12I : lR.oberts et all I2Q11I : 
iPiran et al.l l2Q12f ) which could potentially be observed 
by current and upcoming wide-field survey telescopes 
like PTE (Rau 2009), Pan-STARRS ([K aiser 2004), and 
LSST (|LSST Science Collaborations et al.ll2009[ ). 

There have been nu merous studies o f prim ordial binary 
NS mergers (see e.g. Fa ber fc Rasiol (|2012l )). which will 
have essentially zero orbital eccentricity when they enter 
the frequency band of ground-based GW detectors. How- 
ever, binaries may also form via n-body interactions in 
dense stellar regions and some fraction of them will have 
sizable eccentricity at merger. A likely environment to 
find such binaries is gl obular clusters (GCs) that have un- 
dergone c ore collapse (iFabian et al .11 19 751 : iGrindlav et al.l 
I2006D . In lLee et all (j2010D . it was argued that the rate of 
tidal capture and collision of two NSs in GCs (using M15 
as a prototype) peaked around 2: 0.7 at values of ~ 50 
yr~^ Gpc~^ (falling to ~ 30 yr~^ Gpc~^ by z = 0) and 
was consistent with the observed SGRB rate. However, 
this does not take into account natal kicks. A recent 
simulation of Ml 5 that assumed a modest NS retention 
frac tion of 5% found ^ 1/4 fewer NSs in the central 0.2 
pc (jMurphv et al.l 120111) compared to an earlier study 
that ignored kicks (jPuh et al.lll997L I2003D . implying the 
above rate (scaling as the number density squared) could 
be overestimated by an order of magnitude. On the other 
hand, observations suggest that the NS retention fraction 
in some GCs can be as larg e as 20 % (Pfahl et al . 2002), 
and the model of iLee et"al] ()2010f ) did not take into ac- 
count other channels that could lead to binary merger 
within a Hubble time, suc h as Kozai resonance in a triple 
system (jThompsonl l20 1 if ) . 



The above discussion focused on GC environments, 
and similar i nteractions in galacti c nuclei would also add 
to the rates dO 'Learv e t"al] 120091 : iKocsis fc Levinll20Tl 
lAntonini fc Per ets 2012). Still, it is far from certain that 
high eccentricity mergers occur frequently enough to ex- 
pect observation with the upcoming generation of GW 
detectors. However, it is also not implausible that they 
do, and as eccentric NS mergers may also produce distin- 
guishable EM emission compared to quasi-circular merg- 
ers, it behooves us to understand both systems from a 
multi -messenger pers p ective . 

In iStephens etHI (j2011f ) and lEast et al.l (j2012bD . 
black- hole-neutron- star (BH-NS) mergers formed 
through dynamical capture were found to exhibit a rich 
variation with impact parameter, in some cases pro- 
du cing sizable disks and amounts of unbound material. 
In iGold et al.l ()2011l ). several eccentric NS-NS mergers 
were studied using a F = 2 equation of state (EOS) 
and shown to exhibit /-mode excitation during close 
encounters. There have also been studies of B H-NS and 
NS-NS collisions with N ewtonian gravity (|Lee et al.l 
I2OIOI : iRosswog et al.l l2012f ) showing similar variation in 
the outcomes. 

In this Letter, we study dynamical capture NS-NS 
mergers for a range of impact parameters using general- 
relativistic hydrodynamics (GRHD). We also consider 
several different NS EOSs because of the uncertainty re- 
garding the correct description of matter above nuclear 
densities. One of the important issues we address for the 
first time is if these mergers can produce hypermassive 
neutron stars (HMNSs). In studies of quasi-circular sys- 
tems it was found that thermal energy from the merger, 
as well as differential rotatio n, could support long- lived 
HMNSs for some EOSs fe.g.. iSekiguchi et al.llMlh and 
that this would be imprinted in the GW signal and re- 
sulting disk properties. HMNSs with longer lifetimes 
can also build up significant magnetic fields which can 
pow er strong EM tran sients during the collapse to a 
BH (jLehner et al.ll2QTTI ). For dynamical capture bina- 
ries, the amount of angular momentum, and likely the 
amount of shock heating, will be strong functions of im- 
pact parameter, suggesting HMNS formation will be as 
well. 
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Another notable feature of dynamical capture NS- 
NS mergers is their potential to produce unbound nu- 
clear material which will decompress and form heavy 
nuclei via the r-process (|Lat timer fc Schramml 119741 : 
iRosswQg et"an 119981 : iLi fc Paczvnskil 0998^ ; subsequent 
radioactive decay could produce observable emission . 
Recent work (jFischer et al.|[2QTQl : lArcones fc Jankall2Qllf ) 
suggests processes like NS-NS mergers may be needed 
to supplement the supernovae r-process yield in ac- 
counting for the observed abundances. Though sim- 
ulations of quasi-circular NS-NS mergers using New- 
tonian or conformally flat gravity have found suitable 
ejecta, they seem to be in tension with fully general- 
relativi stic results which f ind negligible amounts of 
ejecta (jFaber fc Rasiol [2Q12f ). This is arguably because 
of strong- field GR effects, such as BH formation and the 
existence of innermost stable orbits. As we show, dy- 
namical capture mergers are more promising sources of 
ejecta, presumably as the stars are less bound when dis- 
ruption occurs. 

In the remainder of this Letter, we outline our methods 
for simulating NS-NS mergers with GRHD, discuss the 
merger dynamics for a range of impact parameters and 
three different EOSs, and comment on potential GW and 
EM counterparts. We find that, while the GW signals 
from these mergers may be challenging to detect with up- 
coming ground-based detectors, they have the potential 
to source numerous EM transients. Non-merging close 
encounters can induce tidal deformations strong enough 
to crack the NSs' crusts; a merger where the total mass 
is above the maximum mass of a single NS can either 
promptly collapse to a BH or produce a hot, rapidly ro- 
tating HMNS, where the latter outcome tends to have 
more massive disks and ejected material. 

2. NUMERICAL APPROACH 

We numerically solve the Einstein equations, dis- 
cretized with finite differences, in the generalized har- 
monic formulation. The hydrodynamics are evolved in 
a conservative formulation using high-resolu tion shock- 
capturi ng techniques. Details are given in lEast et aTl 
([20123 ). 

We use adaptive mesh refinement with up to seven lev- 
els that are dynamically adjusted according to trunca- 
tion error (TE) estimates. To measure convergence and 
TE, we perform a select number of simulations at three 
different resolutions. The low, medium, and high reso- 
lutions, respectively, have base levels covered by 129^, 
201^, and 257^ points (with the maximum TE threshold 
adjusted accordingly), and approximately 64, 100, and 
128 points across the diameter of the NSs on the finest 
level at the initial time (for the HB EOS). In Figure [U 
we show an example of convergence of NS trajectories as 
well as the constraints of the field equations. All simu- 
lations are performed at medium resolution and results 
quoted below are from this resolution, with Richardson 
extrapolated values given in parenthesis (indicating the 
quantity's TE) where multiple resolution data are avail- 
able. 

We use the same gau g e, slop e limiters, and flux 
methods as in lEast et all (j2012b[ ). We use the piece- 
wise polytropic EOS m odels labeled 2H, HB, and B 
from I Read et"al] ()2009l ) and include a thermal compo- 
nent Pth = (Eth - l)ethp with Fth = 1.5. These EOSs 
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Fig. 1. — Top: L^-norm of the constraints (Ca = Ha — ^Xa) in 
the lOOM X lOOM region around the center of mass in the equa- 
torial plane for the rp = 7.5, HB case, with the three resolutions 
scaled assuming second-order convergence. Bottom: the difference 
in the NS center of mass as a function of time from the different 
resolution runs for the rp = 10, HB case, scaled assuming second- 
order convergence. 

were designed to span the range of possible EOSs. The 
2H, HB, and B EOSs, respectively, give NSs with com- 
pactions M^s/Rns of 0.13, 0.17, and 0.18 for Mns = 1.35 
Mq and maximum masses of 2.83, 2.12, and 2.0 Mq 
(unless otherwise stated we use geometric units with 
G = c = l). 

We construct initial data by solving the constraint 
equations in the conformal th in-sandwich formulation as 
described in lEast et al.l ()2012f ). We begin the two NSs at 
a separation of d = 50M, where M = 2.7 Mq is the to- 
tal mass of the system (hence d = 200 km), and consider 
various initial velocities which we label by r^, the peri- 
astron distances of parabolic Newtonian orbits with the 
same velocities (which will be different from the actual 
periastron distance of the simulated binaries). We per- 
formed the majority of the simulations using the middle 
compaction HB EOS but ran select impact parameters 
using all three EOSs. For all simulations except one we 
used NSs that both have a mass of 1.35 M©; the other 
case has a mass ratio of q = 0.8 and a total mass of 
2.88 Mq. We leave a more detailed study of mass-ratio 
dependence to future work. 
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Fig. 2. — Real (solid black) and imaginary (dotted red) compo- 
nents of the Newman-Penrose scalar ^4 on the axis orthogonal to 
the orbit measured at r = lOOM for the rp = 5, 7.5, 8.75, and 10 
cases with HB EOS. 

3. RESULTS AND DISCUSSION 
3.1. Effect of Impact Parameter 

Using the HB EOS, we consider a range of impact 
parameters from r^/M = 2.5 (we henceforth quote 
in units of M) to = 20 (i.e., 10 to 80 km). This 
is weh within the range to form a bound system as 
iPeters fc Mathewsl (|1963f ) indicate that for equal masses, 
GW capture occurs for < l.^/w^^^ where w is the ve- 
locity at infinity. (Tidal energy loss (|Press fc Teukolskvl 
Il977[ ). because of the relative scalings with distance, is 
subdominant in determining capture.) Binaries that ap- 
proach with small impact parameters (r^ = 2.5 and 5) 
promptly merge and collapse to a BH. For = 2.5 the 
mass and dimensionless spin of the final BH is Mbh/^ = 
0.998(0.995) and a = 0.537(0.538) while the energy and 
angular momentum in GWs is EQ^^/M = 3.7(4.0) x 10~^ 
and Jgw/M^ = 2.60(2.75) x 10"^ For rp = 5, 
Mbu/M = 0.985, a = 0.719, Eqw/M = 1.06 x 10"^ 
and Jgw/M'^ = 6.74 x 10"^; Figure [2] shows the corre- 
sponding GW signals. For both these cases, the amount 
of material leftover after merger is < 10~^ the original 
rest mass of the NSs. The dearth of matter post-merger, 
and the fact that most of the power of the GW signal is 
at a relatively high frequency (~ 5 kHz), makes these sce- 
narios less promising sources of observable EM or gravi- 
tational radiation. 

Binaries with larger impact parameters (r^ = 10, 15, 
and 20) result in non-merging close encounters followed 
by long elliptic orbits which we did not follow in their en- 
tirety due to limited computational resources. The close 
encounters result in pulses of GW radiation and excite 
/-mode oscillations in the stars, which are also evident 
in the GW signal (see Figures [2] and [3]). This /-m ode 
excitation was studied in detail in lGold et al. (2011' ) and 
also found in sim i lar BH -NS encounters ptephens et al.l 
120111 : lEast et al.l l2012b[ ). The induced tidal elliptic- 
ity in the = 10 case is greater than the 6R/R ^ 
0.1 value required to induc e a strain of i^strain — 0. 1 
and shatter the NS crust ([Horowitz fc Kadaul (|2009f ): 
though we are not modeling the crust here). The en- 



ergy and angular momentum radiated in the = 10 
close encounter is Eqw/M = 1.472(1.474) x 10"^ and 
Jgw/M^ = 3.545(3.546) x 10"^; for rp = 15 and 20 
(£;gw/M, Jgw/M2) = (1.64 X 10-^^,8.69 x 10"^) and 
(3.8 X 10-^2.8 X 10-^), respectively. Taking this as or- 
bital energy and angular momentum loss gives a Newto- 
nian estimate for the time to the next close encounter of 
65 ms for the = 10 case. For the next largest impact 
parameter simulated, rp = 15, the tidal deformation is 
negligible, and the estimated time to the next close en- 
counter is 1.8 s. This suggests precursor EM transients 
associated with crust shattering for these systems could 
be produced of order hundreds of milhseconds, but prob- 
ably not more than a few seconds, before merger. 

For the intermediate cases (r^ = 7.5 and 8.75), the 
stars come into contact and form a single object. For 
rp = 7.5, this happens at the first close encounter; 
for rp = 8.75, the stars briefiy fiy apart before merg- 
ing. When the stars come into contact they undergo 
shock heating and develop feat ures similar to t h e Kelv in- 
Helmholtz vortices observed in lRosswog et al.l ()2012l : see 
Figure [3]) . Though the total mass is above the maxi- 
mum for a cold, static star with this EOS, the stars are 
highly spun-up and have a significant thermal component 
(22%-25%) to their internal energy (see Tabled]). In the 
vicinity of these objects the (density- weighted) average 
thermal specific energy is eth ^ 10-20 MeV/rrin where 
rrin is the neutron mass. 

These HMNSs produce quasi-periodic GWs with fre- 
quency ^3.2 kHz (see Figure [2|). Though these hyper- 
massive configurations survive the duration of our simu- 
lations (^ 13 ms), they presumably will eventually col- 
lapse to form BHs. In Table [TJ we indicate the rate of 
energy and angular momentum loss to GWs at the end of 
the simulation. From this one can roughly estimate the 
time it will take for the HMNS to radiate its remain- 
ing angular momentum to GWs assuming a constant 
dJcw/dt (e.g., for the rp = 7.5, HB case, it will take 
~ 70 ms). However, magnetohydrodynamic effects, such 
as the magnetorational instability, as well as cooling by 
neutrino emission, none of which we take into account, 
will also be important in determining when these stars 
collapse. Table [T] also lists data from the unequal mass 
ratio {q = 0.8), rp = 7.5 case, which shows qualitatively 
similar behavior to the equal mass case. 

3.2. Effect of Equation of State 

In addition to the HB EOS, we also simulated an inter- 
mediate impact parameter rp = 7.5 using the B (softer) 
and 2H (stiffer) EOSs. For the B EOS, a BH forms 
soon after merger with Mbyl/M = 0.988 and a = 0.766; 
the total radiated energy and angular momentum are 
Egw/M = 3.55 X 10-2 and Jqw/M^ = 0.239, respec- 
tively. For the 2H EOS, the total mass is below the 
maximum for a stable cold NS and the stars fiy apart 
after the first collision before eventually settling down to 
a single massive star. 

We also performed simulations with the B and 2H 
EOSs and rp = 10. With the softer B EOS, the NSs 
undergo a close encounter that is qualitatively similar to 
the HB EOS. However, because of the greater compact- 
ness of the NS, the amplitude of the resulting GW pulse 
is larger with 19% and 13% more energy and angular 
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Fig. 3. — Snapshots of thermal specific energy (top left panel) with a logarithmic color scale from 0.05 to 50 MeV/mn and rest-mass 
density (other five panels) from 10^ to 10^^ gm cm""^. The top left and top middle, and bottom middle and bottom right panels show 
the equatorial plane, while the other two show a perpendicular plane through the center of mass. The top panels show an HMNS with 
surrounding disk and unbound material from the Vp = 8.75, HB EOS case at t = 13.3 ms. The bottom panels show, from left to right, a 
BH and a surrounding disk for the = 7.5 B EOS case dX t = 10.2 ms; a merger from the Vp = 7.5, HB EOS case at t = 3.2 ms; NSs 
with excited /-mode perturbations post close encounter from the rp = 10, HB EOS case at t = 4.1 ms. The first four panels have the same 
distance scale, where the coordinate radius of the HMNS and BH are ~ 13 and ~ 6 km, respectively. The last two panels share a second 
distance scale; the coordinate separation between the NSs in the last panel is ~ 73 km. 
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TABLE 1 

Properties of Hypermassive NS Cases, Measured at t ^ 13.3ms 



^Global angular momentum. 

^Density- weighted average of the thermal component of the specific energy in units of MeV/mn- 
^Fraction of internal (Eulerian) energy that is thermal. 
•^Rest mass that is unbound in percent of Mq. 

®The total energy emitted in GWs through the r = lOOM surface. 

^The total angular momentum emitted in GWs. 

^Average GW fiux of energy. 

^Average GW fiux of angular momentum. 



momentum, respectively. Due to the large eccentricity, 
this will have a significant effect on the subsequent orbit 
of the NSs. The estimate of the time to the next close 
encounter is 50 ms with the B compared to 65 ms for the 
HB EOS. Binaries with impact parameters in this range 
may undergo multiple close encounters before merging, 
with the time between the GW bursts a sensitive function 
of (in this example) the EOS. For the 2H EOS, the stars 
have larger radii and graze during the close encounter, 
merging to a massive star soon after. 

3.3. Possible Post-merger Transients 

Intermediate impact parameters (r^ = 7.5 for the vari- 
ous EOSs and mass ratios, rp = 8.75 HB, and rp = 10 2H 
EOS) form HMNSs with non-negligible accretion disks, 
thought to be necessary for an SGRB progenitor, and un- 
bound material which could potentially power other EM 



transients (as would presumably a subset of larger im- 
pact parameter systems were we to follow them through 
merger). Figure |4] shows the amount of matter, total 
and unbound (fiuid cells with outward radial velocity and 
four- velocity time component Ut < —1; see also Table [1]), 
outside a given radius from the center of mass and the 
velocity distribution of the unbound matter. The var- 
ious cases have 0. 005-0. 07M0 unbound material, and 
roughly 2-3 times more in a disk. As expected, cases 
with less compact NSs have more unbound material com- 
pared to more compact cases. Larger impact parameters 
(which have more angular momentum) also have compar- 
atively more unbound material with the most occurring 
in cases where the NSs first come into contact in non- 
merger close encounters (r^ = 8.75 HB and rp = 7.5 
and 10 2H). The unequal mass-ratio merger with q = 0.8 
produces slightly less unbound material than the com- 
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Fig. 4. — Top: the total and unbound rest mass outside a given 
radius from the center of mass for various cases at t 13 ms. 
Bottom: unbound rest mass with asymptotic velocity grouped in 
0.05c bins. The legend applies to the top panel as well. By this 
time, the rp = 7.5, B case has collapsed to a BH, the rp = 7.5, 2H 
case is an NS below the maximum mass for this EOS, while the 
rest are HMNSs. 



parable equal mass merger. In all cases, the eject a is 
mildly relativistic with asymptotic velocity that peaks 
in the range 0.1-0. 3c. This ejecta is presumably neu- 
tron rich and will convert to heavy elements through 
the r-process, the heaviest of which wi ll undergo fission, 
emitting photons (iLi fc Paczvnskil 1199 8: Kulkarni 2005; 
iMetzger et al.ll20"TQ[ ). The arguments from .Metzger et al.i 
()2010f ) estimate the time scale as 

Wk ^ 0.6 d(Mu/3 X lO-'M0)i/2(^/O.2c)-i/2 
with a luminosity, peaking in the optical/near UV, of 
L ^ 4 X 10^^ erg s-\M^/3 x 10'^ Mq)^^^ {v/0.2cy^^ 

normalized here to the approximate values from the 
rp = 8.75 case. However, recent calculations using 
more detailed heavy element opacities suggest that the 
timesca le may be up to a week with emission peaking in 
the IR (|Kasenll2012h . 

This ejecta is also expected to collide with the in- 
terstellar medium producing radio waves that will peak 
on timescales of weeks with brightness (;Nakar fc Piraui 



Imil) 

F(i/obs)«0.4(£kin/2 X 10^1 erg)(no/0.1cm-3)7/8 

(t;/0.2c)"/^(jyobs/GHz)-3/4((i/l00Mpc)-2 mJy 

where no is the density of the surroundi ng; environment 
(we u se Ho ~ 0.1 cm~^ for GC cores; iRosswog et alJ 

Uohs is the observation frequency, d the distance, 
and we have normalized the kinetic energy and velocity 
to the Tp = 8.75, HB simulation. 

4. CONCLUSIONS 

We have performed GRHD simulations modeling dy- 
namical capture NS-NS mergers, giving direct estimates 
of the corresponding GW emission and merger outcome 
varying impact parameter, EOS, and mass ratio. By 
measuring pre-merger tidal deformation and post-merger 
stripped material (bound and ejected), we have also spec- 
ulated on related EM transients. 

Regarding transients that may precede the merger, 
non-merging close encounters can lead to tidal de- 
formations strong enough to crack the NSs' crust 
and t ap into the ^ 10^^ erg stored in elastic en- 
ergy ([Thompson fc DuncanI Il995l ). potentially causing 
fiaring activity from milliseconds up to possibly a few 
seconds before merger. Though a different mechanism 
and time scales, the signature could be similar to reso- 
nance induced crackin g for quasi-circular inspirals pro- 
posed in iTsang et all |20T2). The cracking of the NS 
crust is one possi ble explanation fo r SGRB precursors 
observed by Swift (|Troia et al.ll2010r ). 

We find that dynamical capture mergers can result in 
prompt BH formation or the formation of an HMNS de- 
pending on impact parameter and EOS. The HMNSs will 
be long lived due to their rapid rotation and thermal en- 
ergy, giving them the potential to seed large magnetic 
fields and source intense transients during collapse. 

In contrast to what was found in general-relativistic 
studies of quasi-circular NS-NS mergers, we find that 
dynamical capture mergers can result in massive disks 
even for equal mass binaries, and can result in up to a 
few percent of a solar mass in ejecta. This mildly rela- 
tivistic ejecta can produce potentially observable optical 
and radio transients. The amount of ejecta found here 
is simila r to the 0.009-0.06 M c?) found with Newtonian 
gravity (jRosswog et al.l [20121 ). though not for compara- 
ble impact parameters (r^ < 5). However, what is qual- 
itatively consistent with the Newtonian setups is that 
we observe the largest amounts of unbound material for 
grazing collisions. 

Regarding GW detectability, the high frequency of 
the merger-ringdown or quasi-periodic signals from the 
HMNS win be difficult to observe with AdLIGO. Indi- 
vidual bursts from close encounters would also not be 
detectable except for very nearby events. For exam- 
ple, an Tp = 10, HB EOS merger at d = 100 Mpc has 
sky-averaged S/N for AdLIGO of ~ 0.9. This implies 
that if dynamical capture NS-NS mergers constitute a 
fraction of SGRB progenitors, a further subset of these 
will not have a detectable GW counterpart. GW sig- 
nals from larger Vp binaries undergoing numerous close 
encounters would have larger S/N, and the timing be- 
tween bursts will be a sensitive function of the orbital 
energy, containing information about the EOS, for ex- 
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ample. We d efer a detailed stu dy of GW detect ability to 
future work (jEast et alJl2Q12al ). 

We thank Chris Thompson for useful discussions, as 
well as participants of the KITP "Rattle and Shine" con- 
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